The B c mass up to order aj 

Nora Brambilla and Antonio Vairo 

Theory Division CERN, 1211 Geneva 23, Switzerland 

We evaluate in perturbative QCD, up to order of, the mass of the B c . 
We use the so-called lS-mass in order to improve the convergence of the per- 
turbative series. Our result is E(B c ) pert = 6326lg 9 MeV. Non-perturbative 
effects are discussed. A comparison with potential models seems to be consis- 
tent with non-perturbative contributions of the order —(40 -r- 100) MeV. 
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I. INTRODUCTION 

The discovery of the B c meson (the lowest pseudoscalar 1 So state of the be system) has 
been reported in 1998 by the CDF collaboration in the 1.8 TeV pp collisions at the Fermilab 
Tevatron J[]. The mass has been measured to be 6.40 ± 0.39 ± 0.13 GeV. 

The fact that the mass of the quarks of quark-ant iquark systems built up by b and c 
quarks is much larger than the typical binding energy suggests that these systems are non- 
relativistic, i.e. that the heavy-quark velocity v is small. The typical scales of these systems 
are the binding energy ~ mv 2 and the momentum transfer ~ mv ; moreover, because of the 
non-relativistic nature of the system, m 3> mv 3> mv 2 (for the purpouses of this discussion 
m and v can be identified with the mass and the velocity of the lightest component of 
the bound state respectively). Let us call Aq C d the scale at which non-perturbative effects 
become important. 

If Aqcd <mv 2 , then the scale mv can be integrated out order by order in a s at a scale 
mv ^> fi' ^> mv 2 . The system is described up to order by a potential which is entirely 
accessible to perturbative QCD and at the leading order is the Coulomb potential. Non- 
potential effects start at order a^lnfi' [|],[3|. This kind of system is called Coulombic. Non- 
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perturbative effects are of non-potential type. In the particular situation mv 2 ^> Aqcd 
they can be encoded into local condensates [|J. This condition seems to be fulfilled by the 
bottomonium ground state, which has been studied in this way in || . Also the charmonium 
ground state has been analysed as a Coulombic bound state by the same authors. In both 
cases (but with caveats) the non-perturbative corrections a la Voloshin-Leutwyler, i.e. in 
terms of local condensates, have been claimed to be under control || . 

For heavy quarkonium states higher than the ground state the condition Aqcd <mv 2 
is not fulfilled and non-perturbative terms affect the potential. The system is no longer 
Coulombic. Traditionally the energy of these systems has been calculated within QCD- 
inspired confining potential models. A large variety of them exists in the literature and 
they have been on the whole quite successful (cf. for some recent reviews). However, the 
usual criticisms apply. Their connection with the QCD parameters is hidden, the scale at 
which they are defined is not clear, they cannot be systematically improved and they usually 
contain a superposition by hand of perturbative and non-perturbative effects. For this reason 
a lot of effort has been devoted, over the years, to obtaining the relevant potentials from 
QCD by relating them to some Wilson loops expectation values [§,0- Anyway, these have to 
be eventually computed either via lattice simulations or in QCD vacuum models ||. In the 
specific situation mv ^> Aqcd 3> mv 2 the scale mv can still be integrated out perturbatively, 
giving rise to a Coulomb-type potential. Non-perturbative contributions to the potential will 
arise when integrating out the scale Aq C d- This situation has been studied in ||. We will call 
quasi- Coulombic the systems described by the situation mv 3> Aqcd 3> mv 2 , when the non- 
perturbative piece of the potential can be considered small with respect to the Coulombic 
one and treated perturbation. 

The only available theoretical predictions (to our knowledge) of the B c mass resort to 
(confining) potential models or to the lattice. In this work we will carry out the calculation 
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of the perturbative B c mass up to order a*. We will call it E(B c ) pert ^ This calculation will 
be relevant to a QCD determination of the be ground state if this system is Coulombic or at 
least quasi-Coulombic. Moreover, in the way we are doing the calculation, we also assume 
the T(1S) and the J /if) to be Coulombic or at least quasi-Coulombic systems. The question 
if these assumptions correspond to the actual systems cannot be settled at this point. On 
the other hand there is no a priori reason to rule them out. Let us consider, for instance, 
the argument used in || for the J /if). Lattice data show that the static potential clearly 
deviates from a 1/r behaviour for distances larger than 1 GeV -1 (see and references 
therein). Therefore the J /if) is Coulombic or quasi-Coulombic if the characteristic scale of 
the bound state, fi ~ m c v c , is bigger than 1 GeV. If we assume m c ~ (1.6 -v- 2.0) GeV and if 
we fix that scale on the Bohr radius, a, of the J/if), \x = 2/a(/i), then we get fi ~ (1.5 -j- 1.6) 
GeV. Since this scale falls into the energy window between the mass scale m c and 1 GeV, 
these figures are consistent with a Coulombic or quasi-Coulombic picture of the J / if) Q. 

The main problem of the calculation of the perturbative B c mass is the well-known 
bad convergence of the perturbative series when using the pole mass. This is due to a 
renormalon cancellation occurring between the pole mass and the static Coulomb potential 
lQfl. We handle the problem by expressing the c and the b pole mass in the perturbative 



It is in general somehow ambiguous to separate in a physical quantity perturbative from non- 
perturbative contributions. From this point of view the following Eq. (||) may be seen as a definition 
of what we call here perturbative B c mass. Analogous definitions have to be understood for the 
perturbative J /if) and T(IS') masses. 

2 One may wonder if these figures are consistent with the non-relativistic expansion underlying 
NRQCD. Only an actual calculation may decide this, since a break-down of the NRQCD expansion, 
if it occurs, should be manifest in a breakdown of the expansion of the energy levels. In the specific 
situation of the B c , as we will see later on in this paper, the expansion that we get shows a still 
convergent behaviour. 
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expression of the B c mass as half of the perturbative mass of the J /if) (E(J/if>) peT t) and the 
T(IS') (-E'(T(lS'))pert) respectively. This corresponds to using the quark mass in the so-called 
15 scheme introduced in [O. In this way, by expressing E(B c ) pert in terms of quantities 
that are infrared safe at order Aqcd (the 3 Si perturbative masses), the pathologies of the 
perturbative series, due to the renormalon ambiguities affecting the pole mass, are cured. 
We will explicitly show that, in fact, we obtain a better convergence of the perturbative 
expansion and a stable determination of the perturbative mass of the B c , just in the energy 
range that the above discussion on the J /if) suggests to be also the relevant one for the B c . 
Non-perturbative terms are of potential type in the quasi-Coulombic situation and of non- 
potential type in the Coulombic situation. They affect the identification of the perturbative 
masses E(J/ip) pevt , E(T(lS)) pevt and E(B c ) pert with the corresponding physical ones. If we 
aim at obtaining a good estimate of the physical B c mass, it is not important for each of 
these contributions to be individually small, as long as the sum of them in the B c mass is 
small. As we will discuss at the end, a picture with non-perturbative corrections to the B c 
mass of a not too large size ( < 100 MeV) seems to be consistent with the experimental data 
and with the potential models. 

The paper is organized in the following way. In the next section we set up the formalism 
and perform the calculation of the perturbative B c mass. In section 3 we briefly discuss the 
non-perturbative corrections and compare our result with other determinations of the B c 
mass available in the literature. 

II. CALCULATION OF E(B c )pert 

In order to calculate the B c mass in perturbation theory up to order a*, we need to 
consider the following contributions to the potential: the perturbative static potential at 
two loops, the 1/m relativistic corrections at one loop, the spin- independent 1/m 2 relativistic 
corrections at tree level and the 1/m 3 correction to the kinetic energy. We will not consider 
here effects due to a non-zero charm quark mass on the b (and the bb system) of the type 
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discussed in [|12j . We will follow the derivation of the heavy quarkonium mass of Ref . [0 . 

The static potential at two-loops has been calculated in |jTB|. It is useful, in order to 
perform an analytic calculation, to split it as 

V (r) = v {r) + 5v {r), 



where vq is the part that does not contain logarithms 
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where /3 n are the /3-function coefficients^], c 
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The strong coupling constant a s is understood in the MS scheme. At the scale \i we will 
take the value of a s from the three- loop expression with A^- -3 = (300 ± 50) MeV. The l/m 
relativistic corrections at one loop, the l/m 2 tree-level spin- independent terms and the l/m 3 
correction to the kinetic energy are given by fi"4 
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3 0o = HCa/3 - 2/3JV>, /3i = 34C^/3 - 10N f C A /3 - 2N f C F , ... 

4 The result one obtains by choosing Nf = 4 and A— |^ 4 = 230 MeV has been checked to be 
consistent with the central value and the errors of Eq. (IBI). 



m&, m c and m red = m b m c / {m b + m c ) are the 6, the c and the reduced pole mass respectively. 
Then, the Hamiltonian relevant in order to get the B c mass at accuracy is 

'2 



H(B C ) =m b + m c + 
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Up to order a* the ground-state energy is given by (( ) means the average on the ground 
state): 
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where at leading order 

(5v G c 6v ) = -m red ( ^ -- ^-ln(/io/2) + 7 ln 2 (/ia/2) 
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with a(/i) = l/(m re dC'irQ! s (/x)), the Bohr radius of the system. This corresponds to the 
only contribution relevant at order produced by the Hamiltonian ([]]) in second-order 
perturbation theory (G c stands for the Coulombic intermediate states) and can be read off 
from the second reference in [El. The other averages can be easily evaluated by means of the 
standard formulas 



- ln 2 (/ir)\ = - ( ln 2 (/ia/2) + 2(1 - lE ) ln(/ia/2) + (1 - lE ) 2 + ^ 
r / a \ b 
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After an explicit calculation we get from Eq 
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where E Q {jjt) = -m rcd . 

The main problem connected with the perturbative series (|3]) is the bad convergence in 
terms of the heavy-quark pole masses. Let us consider, for instance, \x = 1.6 GeV, m& = 5 
GeV and m c = 1.8 GeV. Then we get E(B c ) pcrt ~ 6149 MeV ~ 6800 - 115 - 183 - 353 
MeV, where the second, third and fourth figures are the corrections of order a^, a 3 and 
respectively. The series turns out to be very badly convergent. This reflects also in a strong 
dependence on the normalization scale \i: at n — 1.2 GeV we would get E(B c ) pert ~ 5860 
MeV, while at fj, = 2.0 GeV we would get E(B c ) pcrt ~ 6279 MeV.Q The origin of this 
behaviour can be understood in the renormalon language. The pole mass is affected by an 
IR renormalon ambiguity that cancels against an IR renormalon ambiguity of order Aqcd 
present in the static potential [jTIJ . The non-convergence of the perturbative series (^) signals 
the fact that large (3q contributions (coming from the static potential renormalon) are not 
summed up and cancelled against the pole masses. A possible solution, in order to avoid 
large perturbative corrections and large cancellations, or, in other words, in order to obtain 
a well-behaved perturbative expansion, is to resort to a different definition of the mass. The 
so-called IS" mass of a heavy quark Q is defined as half of the perturbative contribution of 
the 3 S\ Q — Q mass [jO]]. Unlike the pole mass, the IS" mass, containing, by construction, 



half of the total static energy (2m+V Coul ), is free of ambiguities of order Aqcd- Our strategy 
will be the following. First, we consider the perturbative contribution (up to order af) of 
the 3 Si levels of charmonium and bottomonium: 

E(J/^) pcTt = f(m c ), E(T(lS)) pcrt = f(m b ), 

which are respectively a function of the c and the b pole mass and can be read off from Eq. 
(|) in the equal-mass case, adding to it the spin-spin interaction energy: m(CFa s ) 4 /3. We 



5 The result also depends on the c and b pole masses, which are poorly known. See the following 
discussion. 
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invert these relations in order to obtain the pole masses as a formal perturbative expansion 
depending on the IS" mass. Finally, we insert the expressions m c = J/ip) peT t) and 

m& = f~ l (E(T(lS)) peTt ) in Eq. (|3p. At this point we have the perturbative mass of the B c 
as a function of the J/ip and T(IS') perturbative masses 



E(B c ) pcrt = /(r i (E(j/^) pcrt ),r i (E(T(is)) pert )). 



(4) 



If we identify the perturbative masses E(J/i()) peTt , E(T(lS)) pert with the physical ones, i.e. 
E(J/ip) phys = 3097 MeV and E(T{lS)) phys = 9460 MeV 0, then the expansion (|) depends 
only on the scale \i. 
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FIG. 1. j-E^/E^I an d \E$/E2\ as a function of /x, being E n the order a™ contribution to E(B c ) peit . 

The continuous lines correspond to A— L = 300 MeV, the dashed lines to A— f —~ = 250 and 
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A^ = 350 respectively. E3 vanishes for values of around 1.4 GeV. 
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In Fig. |we show the dependence on /i of the order a\ and contributions to E(B c ) pevt 
respectively. Taking into account that the order contribution vanishes at fi ~ 1.4 GeV, 
the perturbative series seems to be reliable for values of fi bigger than (1.2 1.3) GeV and 
lower than (2.6 2.8) GeV. For instance, E(B c ) pcrt = 6278.5 + 35 + 6.5 + 5.5 MeV at the 
scale fi = 1.6 GeV. This is consistent with: i) the fact that, for values of fi close to or less 
than 1 GeV, the perturbative calculation (and the initial assumption that B c is Coulombic 
or quasi-Coulombic) is expected to break down; ii) the fact that higher values of [i do not 
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correspond to the characteristic scale of the system (this is signalled by the appearance of 
big logarithms in the perturbative expansion); in) the estimate of the scale [i inferred in 
the introduction from the size of the J/ip. More precisely we will take in our analysis^ 1.2 
GeV < n < 2.0 GeV and 250 MeV < A^ =3 < 350 MeV (in terms of a s this corresponds to 
0.26 < a s (2 GeV) < 0.30). In this way we entirely cover the energy range used in ||] in order 
to study the J/^-Q 

6.5 
6 . 45 

E(B c ) peTt 6 _ 4 

6 . 35 



ft 

N =3 

FIG. 2. E(B c ) peTt as a function of \i at A^ = 300 MeV (continuous line). The dashed lines 
refer to A^f 3 = 250 MeV and A^ =3 = 350 MeV respectively. 




By varying /i from 1.2 GeV to 2.0 GeV and A^" from 250 MeV to 350 MeV and by 
calculating the maximum variation of E(B c ) pcTt in the given range of parameters, we get as 
our final result 

E(B c ) peit = 6326+g 9 MeV. (5) 



6 The inclusion of a somewhat higher energy region, which seems to be allowed by Fig. |], would 
not change our final result (||). E.g. taking 1.2 GeV < fi < 2.6 GeV we would get, by keeping the 
same central values as above, E(B c ) peit = 6326^^ MeV. 

7 Actually the range considered in || was 1.36 GeV < \i < 1.76 GeV. 
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The upper limit corresponds to the choice of parameters A^ = 350 MeV, fi — 1.2 GeV, 
while the lower limit to A^~~ = 250 MeV and fi = 2.0 GeV. As a consequence of the now 
obtained good behaviour of the perturbative series in the considered range of parameters, 
our result appears stable with respect to variations of /i (see Fig. |2|). We would like to note 
that the main source of error in Eq. (|5|) comes from the border region 1.2 GeV <fi < 1.4 
GeV at < 350 MeV, where it may become questionable to treat the B c as a Coulombic 

or quasi Coulombic system (see Fig. |l|). 



III. DISCUSSION AND CONCLUSIONS 

We have calculated the perturbative B c mass as defined by Eq. (Q). The problem of the 
bad behaviour of the perturbative series has been overcome by expressing the perturbative 
B c mass in terms of the perturbative J/ ip and T(IS') masses. The series we obtain has a good 
convergent behaviour. This fact is relevant since it shows that the scale hierarchy considered 
in the introduction (m > mv > Aqcd)j which led to the Hamiltonian ([I]), correctly applies 
to the system under consideration.^ In other words, the result we get is consistent with 
the assumption made that the B c system is Coulombic or quasi-Coulombic. Moreover, the 
perturbative series turns out to be weakly sensitive to variations of /i (the renormalization 
scale) and A^ =3 in the range 1.2 GeV < /i < 2.0 GeV and 250 MeV < A^ =3 < 350 MeV. 
The result appears, therefore, reliable from a perturbative point of view. Non-perturbative 
contributions have not been taken into account so far. They affect the identification of the 
perturbative masses E(B c ) pert , E(T(lS)) peTt , E(J/ip) pert , with the corresponding physical 
ones through Eq. (|3|). Let us call these non-perturbative contributions 5E(B C ), SE(T(1S)) 
and 5E(J/ip) respectively. As discussed in the introduction, depending on the actual kine- 
matic situation of the system, they can be of potential or non-potential nature. In the last 



For instance, an analogous analysis carried out on the B s system does not show any sign of 
convergence. 
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case they can be encoded into non-local condensates or into local condensates. There is no 
way to discriminate among these situations, since the size of what would be the energy scale 
mv 2 of the system with respect to Aqcd is unknown. Non-perturbative contributions affect 

the identification of Eq. (|]) with the physical B c mass roughly by an amount ~ ^ 

<)hn ilS )± + SE{B C ). Assuming \6E(J/i/>)\ < 300 MeV, \SE(T(1S))\ < 100 MeV and 



2 

8E(T(1S)) < 5E(B C ) < 5E(J/ip), the identification of our result (|) with the physical B c 
mass may, in principle, be affected by uncertainties, due to the unknown non-perturbative 
contributions, as big as ±200 MeV. However, the different 8E are correlated, so that we ex- 
pect, indeed, smaller uncertainties. If we assume, for instance, 5E(T(1S)) and 5E(J/ip) to 
have the same sign, which seems to be quite reasonable, then the above uncertainty reduces 
to ±100 MeV. Constraining even more the form of SE, by evaluating it from the Voloshin- 
Leutwyler formula (i.e. in terms of local condensates), as given in Ref. ||, we get a negative 
contribution (since the term coming from the J/ip is the dominating one) of less than 100 
MeV. This feature, if preserved also in the other kinematic situations, would confirm, in- 
deed, that the effect of the non-perturbative contributions is not too big and that its effect 
is to lower down the perturbative result given in Eq. (||). More quantitative statements 
are difficult to make, since, differently from the perturbative case discussed in the previous 
section, they appear to be dependent on the choice of the parameters. 

The result we get in Eq. (^j is compatible with the experimental value E(B c ) phys = 
6.40 ± 0.39 ± 0.13 GeV reported in 0. We mention that OPAL reports in ^ 2 candidates 
B c in hadronic Z° decays events, with an estimated mass E(B c ) p ^ ys = 6.32 ±0.06 GeV. Also 
this value compares favourably with ours. Having more precise and established experimental 
data will make it possible to make some more definite statements. In particular, it will be 
possible to give, inside a Coulombic or quasi-Coulombic picture, a precise estimate of the size 
of the non-perturbative effects in the B c mass. In the table, we also report, for comparison, 
some of the other determinations of the B c mass available in the literature. The results 
quoted in (T^-|T9[ rely on potential models (essentially a Coulomb plus a confining potential) 



11 



and are reported without errors. The figure that appears in the table in correspondence of 
Ref. |nj refers to an average of different models performed by those authors. Finally P0| 
reports the result of a very recent lattice calculation. We would like to note that, if one 
assumes that potential models give a B c mass close to reality, then, comparing the potential 
model predictions with Eq. (|), non-perturbative contributions seem to be of the order 
— (40 -T- 100) MeV (consistently with expectations, non-perturbative corrections become as 
smaller as perturbation theory better works, i.e. in correspondence of low values of A^- - 
and high values of fi). Finally, it is interesting to notice that these figures are completely 
consistent with the general discussion on the uncertainties, coming from non-perturbative 
contributions, done above. 



Obtained from 


B c (l l S ) in MeV 


Experiment [FJ] 


6400 ± 390 ± 130 


Eq. (||) (pert, mass) 


6326lg 9 


© 


6264 




6253 




6286 




6255 


© 


6.386(9) (98) (15) 



TABLE I. Mass of the B c meson. The result labelled with Eq. (||) refers to the present work. 
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